Supercurrent Quantization in Narrow Channel SNS Junctions 
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We determine the quasi-particle excitation spectrum in the normal region of a narrow ballistic 
superconductor-normal-metal-superconductor (SNS) Josephson contact. Increasing the effective 
chemical potential in the contact converts the electronic levels into Andreev-levels carrying super- 
current. The opening of these superchannels leads to a supercurrent quantization which exhibits a 
non-universal behavior in general and we discuss its dependence on the junction parameters. 



The Josephson effect a hallmark of superconduc- 
tivity, is of fundamental interest and bears considerable 
potential for applications in superconducting electronics 
Today, the miniaturization of electronic structures 
has reached the regime where the transport proceeds 
via few or even a single conducting channel Q. Us- 
ing gated structures, junctions can be transformed from 
insulating SIS to superconducting SNS states with the 
quasi-particle spectrum evolving from phase-insensitive 
electronic- to phase-sensitive Andreev states carrying 
large supercurrents. The onset of superflow proceeds in 
steps associated with the opening of transverse channels, 
similar to the conductance quantization in a metallic con- 
tact. While in the latter the conductance G ~ I/V is 
quantized in units of 2e'^/h ^]^, it is the maximal super- 
current Ic which is quantized in units of eA/h in the case 
of a short Josephson link (A is the superconducting 
gap in the banks). However, contrary to the universality 
of the quantization in a normal contact, the quantization 
of the critical supercurrent is perfect only in the limit of 
short junctions L <C (Co = hvp/irA is the supercon- 
ducting coherence length) but is non-universal in general 
Indeed, while experiments on superconducting quan- 
tum point contacts do show steps in the critical current 
Ic, these are much less prominent than the corresponding 
steps in the conductance G 0. In this letter, we study 
the opening of superconducting channels in the metallic 
link of a narrow ballistic SNS Josephson contact and de- 
termine the evolution of the quasi-particle spectrum and 
the dependence of the supercurrent quantization on the 
junction parameters. 

While the behavior of macroscopic SNS Josephson 
junctions is well understood the present interest con- 
centrates on gated structures of mesoscopic size. Such 
quantum point contacts are realized in heterostructures 
[^,0 or via manipulations with a scanning tunneling mi- 
croscope 1^. Theoretically, the supercurrent- phase re- 
lation in mesoscopic SNS junctions with a (5-scatterer 
has been analyzed [ic| ] and the phenomenon of super- 
current quantization has been studied in short junctions 

. Non- universal features of supercurrent quantization 
have first been observed by Furusaki et al. — unfor- 
tunately, these numerical results provide limited insight 
into the physical origin and the parametric dependence of 



these effects. Here, we present a detailed discussion of the 
opening of superchannels in mesoscopic SNS junctions us- 
ing quasi-classical and scattering matrix techniques. We 
discuss the non-trivial evolution of the excitation spec- 
trum as the chemical potential drops below the super- 
conducting gap and analyze the transformation of the 
ballistic SNS structure into a SIS tunnel junction. The 
phase-dependence of the quasi-particle spectrum allows 
us to find the supercurrent quantization in short and long 
junctions, where the contribution from the continuous 
spectrum can be ignored. 
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FIG. 1. Narrow channel SNS contact: (a) geometrical 
setup showing the adiabatic joining of the wire to the su- 
perconductor, (b) potential landscape with a flat potential 
in the wire center and smoothly dropping to the band bot- 
tom in the superconducting banks. While electrons and holes 
with small excitation energies e < jj,x form current carrying 
Andreev states, the hole propagation is quenched at large en- 
ergies e > fix and the Andreev levels transform into localized 
electronic states. 

Consider a narrow metallic lead (with few transverse 
channels) connecting two superconducting contacts as 
sketched in Fig. 1(a) [we assume piecewise constant gap 
parameters A(a; < —L/2) ~ Aexp(i(y9L), A(|a;| < L/2) = 
0, and A(L/2 < x) = A exp(i(pR)]. For each transverse 
channel in the metallic wire the quasi-particle spectrum 
El, is determined through the ID Bogoliubov-de Gennes 
equation (we choose states with > 0) 



2m 



fix{x) A{x) 
^*(^) ^+Mx(2:) 







— c- 






v^{x) 







where u^, and denote the electron- and hole-like com- 
ponents of the wave function ^'j,. The spectrum splits 
into continuous and discrete parts and we concentrate 
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on the latter part in the following, < A. 

The effective chemical potential fj.x{x) ~ — ej_{x) 
accounts for the transverse energy £^(2;) of the chan- 
nel, see Fig. 1(b). Within a quasi-classical formulation 
we describe the quasi-particles in terms of their kinetic 
energies K± = h'^k'^/2m — fj,x{x) ± e and assume trans- 
mission and reflection to be ideal (the excitation energies 
£ = E — > are measured with respect to the Fermi 
energy Sp). An electron with energy e < A below the 
gap is reflected back from the superconductor as a hole 
with kinetic energy K- = fj.^ — £, injecting a Cooper- 
pair into the superconducting contact, a process known 
as Andreev reflection jl^ . A second reflection at the op- 
posite NS boundary transforms the hole state back into 
the original electron state, thus producing a phase sen- 
sitive Andreev level carrying the supercurrent across the 
normal metal lead. The hole-part associated with the 
Andreev level can propagate only if its kinetic energy is 
positive, K- > 0, see Fig. 1(b). Otherwise, the hole 
is back-reflected from the normal potential in the junc- 
tion and transformed into an electron at the NS bound- 
ary — the incident electron is effectively reflected back 
as an electron and a phase-insensitive electronic level is 
formed. Hence, the superchannel starts being modifled 
when the chemical potential fix drops below the gap A 
and is quenched completely with all Andreev levels trans- 
formed into electronic ones when fi^ becomes negative. 

Going beyond quasi-classics, the above physics is con- 
veniently described through the scattering matrix for- 
malism 1 12 1^. We deflne scattering states in the 
normal region and characterize them through the en- 
ergy dependent transmission and reflection coefficients 
iexp(ix*) and rexp{ix^) describing the propagation of 
quasi-particles incident from the left through the junc- 
tion. Matching these states with the evanescent modes 
in the superconductors we obtain (within the Andreev 
approximation |l|]: {K+ - K_)/{K+ + K^) <^ I at the 
NS interface) the quantization condition, 

cos(S'+ — S- — a) = r+r- cos (3 + t+t- cos (p, (1) 

where the +(— ) signs refer to the positive and negative 
energies ±£ of the electron(hole)-like quasi-particles and 
S±{e) = X± + ko,±L, with ko^±L = ^/2rn(e^±e)L /Ti the 
phase for free propagation (while the phase S refers to 
the propagation from — L/2toL/2, the scattering phases 
X* and refer to the origin). Andreev scattering at the 
NS boundaries introduces the phase a = 2arccos(e/A) 
ranging from tt at e = to at the gap e = A, as well as 
the phase difference Lp — ^pi^ — Lp-^ between the two super- 
conducting banks. The phase (5 = (x+ — X+)~ (X- 
reduces to /? = for a symmetric barrier in the absence 
of perfect resonances [as follows from the unitarity of the 
scattering matrix; for a symmetric potential shifted by a 
from the center we have (3 = 2(fco.+— fco,-)a]. The secular 
equation (^ involves two main energy dependencies orig- 
inating from the propagation through the normal layer 



(S'lb) and from scattering at the NS boundaries, e.g., due 
to potential steps or barriers. Here, we concentrate on 
the case where the transport through the junction is dom- 
inated by the normal metallic wire — we will comment 
on the effect of resonances introduced by an additional 
scattering at the NS boundaries below. 
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FIG. 2. Discrete energy spectrum [(a) qualitative sketch for 
a flat potential, (b) smooth parabolic potential]: For fj,xo > £ 
(region I) both electrons and holes propagate, forming phase 
(ip) sensitive Andreev levels carrying supercurrent. The dou- 
ble degeneracy of the Andreev states is lifted by a finite phase 
drop ip across the junction as well as a finite refiection in the 
wire (see (c)), the latter becoming relevant upon decreasing 
fixo- As jj,xO drops below e the Andreev levels first transform 
into electronic states (regions II and III) and finally turn into 
boundary states trapped at the NS interfaces when pxo < ~e 
(region V). Within the shaded regions around e = itfixo the 
transmission drops from unity to zero for holes and electrons, 
(c): graphical solution of (1) along the cut A-A in (b). 

A rough understanding of the transformation from 
a metallic to an insulating junction is obtained in the 
quasi-classical approximation using a flat potential, see 
Fig. 1(b): For a large chemical potential fix > s we 
have r± = 0, t± = 1 in (^ and we obtain the Bohr- 
Sommerfeld quantization condition for the (phase sensi- 
tive) Andreev levels 5*+ — S- ~ a ± p — 2mr. Evalu- 
ating this condition for a flat potential (ignoring contri- 
butions to S± = kp^xLy^l ± e/ fix originating from the 
adiabatic joints) we obtain the level scheme shown in 
Fig. 2(a). On the other hand, for — e < fix < £ the 
right hand side of (|l|) vanishes and using S- — 3tt/2 (as- 
suming a hard wall potential) we find the quantization 
condition 25*+ — 2a = 2?T,7r, for the electronic levels, see 
Fig. 2(a) for a qualitative result in a flat potential. Note 
that we have twice as many electronic than Andreev lev- 
els as the latter are doubly degenerate a,t ip = — the 
exact transformation of the Andreev levels in electronic 
ones at e ~ fix requires a more careful analysis account- 
ing for the non-ideal transmission and reflection through 
the normal channel, see below. Finally, as fix drops be- 
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low — e both the electron- and hole-like trajectories are 
quenched. Note the evolution of the r.h.s. of (|l|), going 
from cos Lp at large positive > e, to a small value in 
the intermediate region — e < /i < e, and back to unity 
at large negative ^l^ < —s. This provides us with a first 
rough understanding of the SNS to SIS transformation. 

In a more accurate study of the evolution of the bound 
state spectrum from a SNS to a SIS junction we as- 
sume a smooth potential fj.x{x) with a small curvature 
mfl^ = dxfJ'x, < A, producing a sharp switching be- 
tween transmission and reflection within the energy in- 
terval hfl (a ^-function scatterer Q does not describe a 
pronounced transformation from a SNS to a SIS junc- 
tion). Adiabatic joining of the wire to the supercon- 
ducting banks requires that TOr2^(L/2)^/2 ^ £p and al- 
lows us to make use of the Andreev approximation while 
avoiding the appearance of resonances (this condition 
can be relaxed as the Andreev approximation requires 
mQ'^{L/2)'^/2 > A, while a step in the potential AV" < 
0.9£f produces only weak resonances]. In summary, a 
smooth contact without resonances requires the param- 
eter setting ^SpEi^ < hn < A where = fi^TT"^ /2mL^; 
this condition implies a long junction L > £^nd hence 
a relatively large number n ^ ^Je-p/ A. of trapped lev- 
els. For such a smooth potential the Kemble formula 
is valid and the transmission probabilities take the form 
T±^tl^ 1/{1 + exp[-2^(Ai,(0) ± e)/nn]} %. 

Fig. 2(b) shows the refined results for the SNS 
to SIS transformation using the quadratic poten- 
tial Hx{x) = HxQ + with the parameter 
hn{pxo) = (4/7r)A/£i(eF - Uxo) and ^ixo = ^J■xiO), 
joining the band bottom of the superconductors at 
the two NS boundaries. For this case, the quasi- 
classical dimensionless action takes the form S{E)/h ~ 
{2E/nn)[K^^/l + K-2 -I- ln[|K|(l + \/TT^)] with k2 ^ 
Qhfl/E — TT^h^Cl^ /IGEel, Q » 1 a large parameter 
[S± = S{E = fixo i e); an additional phase tt, which 
cannot be obtained within the quasi-classical scheme, is 
picked up over the energy interval hn as E goes through 
zero]. As for the flat potential, the Andreev levels at 
large chemical potential fixo > £ + hQ (region I) are con- 
verted in steps (regions II — V) to the electronic states 
at negative potential fixO < — £ ~ Upon entering re- 
gion II the product in (0) drops below unity and the 
Andreev levels split even for ip = 0. In region III the hole 
propagation is quenched and electronic levels with an ex- 
ponentially weak phase sensitivity are formed. Entering 
region V the electronic propagation through the wire is 
suppressed and pairs of boundary states are formed near 
the NS interfaces (for the flat potential these bound- 
ary states are lifted to e sa A as the Andreev phase 
a has to vanish). In region IV both the electron- and 
hole-components undergo flnite reflection and the dis- 
tinction between Andreev- and electronic levels is gone; 
this resembles the situation of a SNS junction with a S- 
scattering potential Cn] . More details of this conversion 



between Andreev- and electronic levels will be given else- 
where . Below, we concentrate on the quasi-classical 
region I and study the evolution of the critical supercur- 
rent as the channel is switched on and off. 

The supercurrent / flowing through the junction splits 
into the two contributions from the discrete (Jdis) and the 
continuous (/con) parts of the spectrum. Here, we con- 
centrate on /dis, which dominates the expression for the 
critical supercurrent in the quasi-classical regime I (re- 
quiring hil < iJ,xo is sufficient). The current of the ly-th 
level (including a factor 2 for spin) can be obtained from 
the derivative Ii, = {2e/h)d^ei, = (2e/T) t^t^ sinip, with 
the generalized travelling time T = shi(6S^a)hde [SS—a] 
+hd^[{t+t-) cos ip + (r_|_r_) cos/3]. Within the quasi- 
classical region I, each Andreev level carries a finite su- 
percurrent of amplitude 2e/[r-|_ + t_ + 27i/\/ A^ — e'^], 
where t± = TideS± denote the propagation times for the 
electron- and hole-components; for the smooth quadratic 
potential we find t{E) = f]-i[21n[|K|(l + 71+"^)]. 
For small energies the travel time increases logarith- 
mically t{E) w Vl^^ \a.{iQhn / E) within the interval 
fiQ. < E < QMl and saturates at tq « ln(4(5) as E 
drops below Tifi, a result going beyond the quasi-classical 
approximation. Kt ip — Q, the pairwise degenerate lev- 
els produce equal currents of opposite sign and the sum 
over the discrete spectrum gives no current. Increasing 

the degeneracy is lifted and the resulting miscancella- 
tion leads to a finite supercurrent. Each pair produces a 
monotonuously growing current of the same sign, hence 
the largest current is reached at ip = tt^ . However, at 
ip> = 7r~ the levels become degenerate again and their cur- 
rents cancel pairwise, except for the lowest level which 
remains unpaired and thus carries all the supercurrent 
from the discrete part of the spectrum. The continuous 
part of the spectrum vanishes at </? = tt, however, this is a 
priori not sufficient to guarantee that the critical current 
Ic is the current Iq carried by the lowest level — we have 
to show in addition that the maximum of / = Jdis + Icon 
is reached at p> = tt^ [indeed, we could prove that this 
condition is fulfilled within a regime of the L-^xO plane 
away from (L ^ (<^o/^f)^^^, MkO ~ A)]. In the end, we 
arrive at a particularly simple expression for the critical 
current density in the quasi-classical region I, 

Ic^e/iTo + h/A). (2) 

The travel time tq is constant (il^^ln4Q) at the open- 
ing of the channel, decreases as In AQhil/ ^xo for 
fj-xo > and transforms to the free travel time L/vp^x 
for Uxo > Qhil. As the channel becomes wide open at 
high energies, the critical current saturates to the ex- 
pected value Ic = ev-p/ {L + tt^o)- 

The above discussion dealt with the parameter settings 
^EfEl < Jt-^ < A requiring a junction with L > ^q. Re- 
leasing the condition of small curvature and assuming 
A < hfl, the SNS to SIS transformation is smeared and 
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the region IV occupies all of the interesting crossover 
regime. For a flat potential hi} < s/e^el the situa- 
tion is complicated by the appearance of resonances due 
to reflection from the potential step at the NS bound- 
ary ]l5[ |. The situation simplifies for a very short junc- 
tion with L <C ^Oi where we can again make use of 
(|l|) to produce a simple and universal result with 
J5 = 5+ - S". « and t_K.t+K. ^/T, r- ^ r+ ^ y/R, 
we find that only one level remains trapped in the junc- 
tion at Eo = A[l — Tsin^((/3/2)]^/^. Here, we require a 
width hn > A in order to avoid a strong energy depen- 
dence in the transmission probability T. Determining the 
current Io{ip) from Eq and maximizing, we obtain 

/e = {eA/n) (f - VR), (3) 

in marked difference from the result for the conductance 
quantization G = (2e^/?i)(l — R): a finite reflection < 
R ^ I will affect the supercurrent quantization much 
more strongly than the conductance quantization. 
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FIG. 3. Supercurrent quantization: with increasing width 
d of the normal channel the supercurrent increases in steps 
of e/(ro -I- h/A). The dotted and dashed lines give the ap- 
proximations To « \n{4:QhQ/ E) and tq ~ L/vp at small 
and large energies (parameters: £f = 1 eV, A/ep = 10~^, 
L/^o ~ 1, 10; with the curvature hQ/A < 5, 0.5 no smearing 
is visible at the supercurrent onset). 

Finally, we discuss the supercurrent quantization 
'steps' appearing as the gate potential is decreased to 
open the conducting channel. We concentrate on the 
quasi-classical regime I, assuming a parabolic potential 
fix{x) in the junction which matches the band bottom 
of the superconductors at the NS boundary and ignore 
a possible change in the effective mass. The quantized 
transverse energy of a channel of width d is given by 
e±.i « Tt^'tt'^P /2mcP; these levels match up with the 
Fermi energy when d = dk = kn/k-p. As we open 
the fc-th channel the other open channels have already 
dropped by fJ.xO-i,k = £f - £±-i,k = £f(1 - I'^/k'^), where 
£_L;i,/c = Ti^TT^P /2md1.^ e.g., the first channel is wide open 
when the second channel appears, iJixO;i,2 = (3/4)£f. In- 
creasing the channel width d, the first channel opens 
(i.e., = eF(l — T^'^/k^d'^) turns positive) as we reach 



di = n/kf, the critical current increases sharply Ic ~ 
eil/ ln[4/(A:pd^/7r^ — 1)] (the logarithmic singularity is 
cutoff at Uil) and saturates at Ic ~ evp/{L + tt^q), see 
Fig. 3. Here, we have ignored the smearing near the on- 
set within the range hfl due to a finite reflection — while 
the interesting evolution of the quasi-particle spectrum 
is washed out as Tifi increases beyond the gap A, see Fig. 
2, the (smooth) steps in the onset of the supercurrent 
are much more robust. On the other hand, the absence 
of sharp steps in the critical current onset is an intrinsic 
feature of the superconducting junction. For the short 
junction with L ^ ^ the sharpness of the steps in Ic is 
dictated by the reflection probability R of the junction 
and thus is more similar to the steps in the conductance 
G. However, with /c cx 1 — Vr, the steps in Ic are always 
smoother than those in the conductance G oc {1 — R). 

In the end, universal supercurrent quantization first 
seems to require short junctions, but the gate needed to 
switch the channels will produce backscattering and spoil 
the quantization. While going over to longer contacts 
helps to produce sharp conductance steps, the onset of 
supercurrent remains smooth due to the long travelling 
time for the Andreev states. 
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